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ABSTRACT 

Context. Since 2004 we have been carrying out a pulsar survey of the Cygnus region with the Westerbork Synthesis Radio Telescope (WSRT) 
at a frequency of 328 MHz. The survey pioneered a novel interferometric observing mode, termed 8gr8 (eight-grate), whereby multiple simul- 
taneous digital beams provide high sensitivity over a large field of view. 

Aims. Since the Cygnus region is known to contain OB associations, it is likely that pulsars are formed here. Simulations have shown that this 
survey could detect 70 pulsars, which would increase our understanding of the radio pulsar population in this region. We also aim to expand 
the known population of intermittent and rotating radio transient (RRAT)-like pulsars. 

Methods. In this paper we describe our methods of observation, processing and data analysis, and we present the first results. Our observing 
method exploits the way a regularly spaced, linear array of telescopes yields a corresponding regularly spaced series of so-called "grating" 
beams on the sky. By simultaneously forming a modest number (eight) of offset digital beams, we can utilize the entire field of view of each 
WSRT dish, but retain the coherently summed sensitivity of the entire array. For the processing we performed a large number of trial combina- 
tions of period and dispersion measure (DM) using a computer cluster 

Results. In the first processing cycle of the WSRT 8gr8 Cygnus Survey, we have discovered three radio pulsars, with spin periods of 1.657, 
1.099 and 0.445 seconds. These pulsars have been observed on a regular basis since their discovery, both in a special follow-up programme as 
well as in the regular timing programme. The timing solutions are presented in this paper We also discuss this survey method in the context of 
the SKA and its pathfinders. 

Conclusions. We have found three new pulsars using the WSRT. Reprocessing and further analysis of the data will reveal dimmer pulsars, and 
RRAT-like or intermittent pulsars. 

Key words, stars: neutron - pulsars: general 



^ 1. Introduction 



In the last decade several pulsar surveys have been conducted 
using various telescopes. Many of those surveys have been ex- 
tremely succesful, and altogether they have more than doubled 
the number of radio pulsars known in our galaxy; they have also 
found pulsars in some other nearby galaxies (e.g. Manchester 
let al.|2001[ [Crawford et al.|2001[|CoFdes et al.|,2006, .Hesselsi 



et al. 20081. Ongoing surveys are continuously contributing 



to the understanding of the characteristics of the whole pul- 
sar population such as biith rates, velocity distribution, lumi- 
nosity function, beaming fraction and magnetic field evolution. 
Studying large numbers of sources also provides insight into 
pulsar formation and evolution, and tells us about the local 
density and pressure environments in the interstellar medium 



and the magnetic-field structure and electron distribution of the 
Galaxy itself. Furthermore, unique objects may be discovered, 
opening up new possibilities for pulsar science; e.g. the dou- 
ble pulsar ( Burgay et al.||200? Lyne et aL]|2004 l, or possible 
PSR-BH systems. 

In 2004 we began a special survey using the WSRT in The 
Netherlands. We developed a new beamforming technique that 
makes optimal use of the array configuration and multiband 
IF (intermediate frequency) system of the telescope for sensi- 
tive, wide-field pulsar surveys. We combine the field-of-view 
advantage of the small single dishes with the high sensitivity 
and angular resolution of a tied-array beam. This allows for a 
quick and accurate determination of the positions of all new 
pulsars. This survey covered a region containing the Cygnus 



superbubble and many OB associations (e.g. Uyamker et al. 
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|2001| l. In such regions of high star formation rate, it is rea- 
sonable to expect the discovery of significant numbers of, in 
particular, young pulsars. 

In this paper we describe the setup and characteristics of 
the 8gr8 survey, and present the first results: the discovery of 
three pulsars. 

2. Observations and data analysis 

2.1. 8gr8 survey: general Introduction 

Radio telescopes that are designed to be used as interferom- 
eters are, in general, not useful for large area pulsar surveys. 
This is because they typically consist of several small dishes 
that are later combined, for imaging purposes, using a correla- 
tor. In order to obtain the sensitivity of the full collecting area 
of a synthesis telescope for pulsar observations it is necessary 
to combine the dishes "in phase", which means that the geo- 
metric and instrumental phase terms are corrected before the 
telescope signals are combined. These corrections are typically 
made for only one direction on the sky resulting in a "tied- 
array" beam with a width inversely proportional to the largest 
separation of the dishes in the array. In a typical WSRT obser- 
vation of a known pulsar, all 14 of the 25 m diameter dishes 
are coherently combined to give a collecting area equivalent to 
a single dish of 94 m diameter. However, the size of the resul- 
tant beam on the sky corresponds to a dish with a diameter of 
2.7 km, the largest separation in the array. A simple calculation 
shows that the ratio of the field-of-view of a single 94 m dish to 
that resulting from the 2.7 km baseline of the WSRT is ~ 1000. 
A tied-arTay pulsar survey with the WSRT would therefore take 
~ 1000 times longer than a similar sur-vey using a single 94 m 
dish. 

We have developed a new method that allows us to over- 
come this limitation by making use of the fact that the first 
12 WSRT telescopes, RTO through RTB, can be placed at 
equal separations of 144 m. The linear (East- West) nature of the 
WSRT array means that when the telescopes are added in phase 
to form a tied-array, the instantaneous result is an elliptically 
shaped fan beam with a width (minor axis) inversely propor- 
tional to the separation of the furthest dishes, and a length (ma- 
jor axis) equivalent to the size of the primary beam. The equally 
spaced dishes produce a grating response on the sky, with par- 
allel fan beams spaced at intervals of c/ {B v) radians, where B 
is the projected baseline between the dishes and v is the ob- 
serving frequency (see Fig.|2^). We will refer to the collective 
fan beams in such a grating response as a grating group from 
now on. While these fan beams provide an increased field-of- 
view we can still do better by utilising the fact that the WSRT 
has eight independent signal chains allowing us to form eight 
simultaneous grating groups on the sky (Fig.|2^). Each of these 
grating groups is tuned to the same sky frequency but has dif- 
ferent delay- and phase-tracking centres, offset by c/(B v)/8, al- 
lowing us to almost fully tile the entire primary beam. Although 
the eight instantaneous grating groups do not fully sample the 
primary beam (12 groups would be required for that), the ro- 
tation of the Earth during the ~ 7000 second observations lead 
to smearing of the point spread function, such that the varia- 
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Fig. 1. Combined detection beamplot of the original and con- 
firmation SgrS observations in which PSR J 1937 -1-2950 was de- 
tected. The closed circles represent the subbeams in which the 
pulsar was detected in the first observation, and the open circles 
the detections in the confirmation observation. The sizes of the 
circles correspond to the S/N of the detections. The grating ef- 
fect as described in the text results in multiple detection-regions 
for each observation. Using the two observations at different 
hour angles, centred at the same position, allows for discrimi- 
nating between the regions. The cross indicates the best-fit po- 
sition of PSR J 1937 -1-2950 from the WSRT timing programme, 
see Table[T] 

tions in total integrated sensitivity over the beam are small (see 
Fig.|2]3). The signals from the eight beams are later combined 
(using the appropriate linear combinations) in software to form 
elliptical subbeams which have dimensions equal to the beam 
size that would result from an equivalent synthesis observation. 
These subbeams are spaced to have them overlap at the half- 
power point of both axes. They tesselate out the entire primary 
beam and have a typical sensitivity that can be compared to that 
of a 74 m dish. 

It is important at this stage to note some consequences of 
this method. The size of the primary beam, c/(vD) where D 
is the diameter of the dish, remains constant as a function of 
hour angle. However the separation and width, 1 IB v/c, of the 
fan beams depends on the length of the projected baseline, and 
thus is hour angle dependent. This means that the number of 
fan beams from a particular grating group which fall inside the 
primary beam changes from a minimum of one, when the pro- 
jected baseline is equal to the dish diameter, up to a maximum 
of five at the zenith. 

During an observation, the central fan beams track loca- 
tions close to the primary beam phase centre which means that 
they "rotate" over source locations within the primary beam. 
When observing at large hour angles, the large width of the fan 
beams, combined with our approximately two hour observation 
duration, means that the sources remain predominantly in just 
one grating group. However at observations near the zenith, a 
source can pass through many grating groups. This is compen- 



sated in off-line processing, see Sect. 2.3 As each of the fan 
beams that make up a grating group has almost the same gain, 
any source which appears in one fan beam will also be present 
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Fig. 2. a (left): The instantaneous grating fan-beam response of the WSRT telescopes RTO through RTB, added in phase, when 
all have a relative baseline of 144 m. A cross-cut is displayed through the narrow dimension of the fan-beam. The overall taper 
corresponding to the telescope primary beam has not been included. One of the nearest neighbouring instantaneous grating groups 
is illustrated with the dashed line. A total of eight grating fan-beams provides complete coverage of the primary beam. Although 
the instantaneous fan-beams only overlap at about the 10% level, any finite integration time will broaden the effective beams, 
yielding a more nearly uniform combined sensitivity, b (right): The relative sensitivity across the primary beam field-of-view 
after a 2 hour integration with a total of eight grating fan-beams. The largest oscillations are within a few beamwidths of the 
overall field centre. These damp out at larger distances to a constant level. 



in the other fan beams from the same grating group inside the 
primary beam. An example of this is shown in Fig.[T] This 
means that one 8gr8 observation, with observing parameters as 
described below, does not allow us to distinguish between the 
fan beams from the same grating group when determining the 
position of a new source. 

This observing method has some useful advantages over 
single-telescope, single-receiver observations. While the width 
of the fan beams, or the number of fans from a grating group 
in the primary beam limit the accuracy to which the position 
of a new source can be determined, it is still better than using 
an equivalent observation with a single dish of the same size. A 
corollary of this is that any real source must be detected in more 
than one subbeam. In contrast, any source that appears with 
similar significance in a large number of subbeams is likely to 
be interference. This is because any near-field source should 
appear with equal sensitivity in all 8 of the bands. The almost 
two hour long observation time used for this survey means 
that, depending on hour angle, some positional discrimination 
is possible (see Fig.[T]|. To provide the ability for both direct 
confirmation of candidate pulsars and to improve positons for 
those candidates, we made a second observation of all point- 
ings at different hour angles, leading to fan beams with differ- 
ent position angles within the primary beam. When combined 
with the first pointing this allows for further restriction of the 
likely position of the source. An example is shown in Fig.[T] 

2.2. Observations 

The region we have surveyed extends along the Galactic plane 
40° < / < 100°, -0^25 < b < 7°25, covering a total of 450 



square degrees. The survey was carried out at a central fre- 
quency of 328 MHz with a bandwidth of 10 MHz. The field- 
of-view of the primary beam of the 25 m dishes at the WSRT at 
this frequency is 5.2 square degrees and, allowing for a small 
degree of overlap between the beams, we were able to span this 
region using 72 pointings spread over 12 days. We observed 
for 6872 s (2^^ samples of 819.2yus) and to be able to reach 
the high DMs that might be expected in the Cygnus region, 
we used 512 channels across the 10 MHz bandwidth. Data rate 
limitations in the Pulsar Machine (PuMa) observing hardware 
( Voute et al.|2002| l meant that we had to reduce the sampling 
interval to a modest 819.2jt/s, meaning that our sensitivity to 
the fastest rotating pulsars was reduced. However, as we were 
expecting, and aiming to detect, predominantly young pulsars 
in this region of the sky, this compromise was considered to be 
optimal. Motivated by the possibility that nearby pulsars might 
be missed in our survey due to scintillation, we observed each 
of the 72 pointings twice to increase the chance of detecting 
the source in a scintillation maximum. As mentioned before, 
this second data set also enabled instant confirmation of new 
pulsars. 

The maximum sensitivity of the survey is shown in Fig.|3] 
using the survey parameters as described above, and represents 
an eight sigma detection. We consider three different Galactic 
latitudes, b - 0, 3, 6 degrees, to indicate the increase in sensi- 
tivity as Tsky decreases when we move away from the Galactic 
plane. We calculate the sensitivity for two representative DMs: 
10 and lOOcm^^pc. The calculation also includes the contri- 
bution to pulse broadening from scattering, as estimated using 
the DM-scatter broadening relation of ,Bhat et al.| ( j2004^ ). These 
curves represent the maximum sensitivity at the pointing cen- 
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( [Manchester et al.|19"96i [Edwards et al 



2OOTI1. A modified form 



Fig. 3. Sensitivity curves, calculated using Eq. A1.22 of 
[Lorimer & Kramer[ ( |2005| l, are shown for DMs of 10 and 100 
and Galactic latitudes b - 0,3,6 degrees, at the observing fre- 
quency of 328 MHz. The curves correspond to the minimum 
flux for detections of 8 cr. 

tre and should be convolved with the curve shown in Fig.|2|3 to 
give the sensitivity across the field-of-view. 

Where possible, observations were scheduled at large hour 
angles so as to maximise the width of the fan beams and thus 
minimise the number of sub-beams that needed to be formed 
to tile out the primary beam. Despite this, the number of sub- 
beams ranged from a minimum of 400 to a maximum of 2300 
for the different pointings, resulting in a total of approximately 
100 000 subbeams across all 72 pointings. As each pointing 
was observed twice, and each subbeam has to be formed and 
then searched in the usual way individually, this results in a 
very large computational load. 

The initial observations were completed in 2004 and the 
second set was finished in 2005. The massive computational 
load of searching more than 200 000 subbeams meant that 
searching of the data could not begin until the arrival of the 
PuMa II cluster ( [Karuppusamy et al.| 2008 ) in early 2005. The 
need for that cluster to be used for other observing programmes 
has meant that to date only a relatively low sensitivity pass 
through the data has been possible. In order to improve process- 
ing speed, the data were initially processed without any time 
domain interference removal nor correction for any baseline 
variations. The latter occurs when bright interference causes a 
readjustment in the level settings for packing the data into two 
bits. It results in steps in the data on timescales of tens of sec- 
onds which reduce our sensitivity to pulsars with periods longer 
than about a second. We have developed a new processing pro- 
cedure which is able to correct for these steps and also is more 
robust to time domain interference. The significantly improved 
analysis, including single pulse searches, is currently ongoing 
and will be presented in a forthcoming paper 



of the tree dedispersion algorithm ( Taylor] 1974 1 was applied to 
each of the data sets from the 8 beams. The total bandwidth of 
each beam was divided into smaller frequency subbands within 
which the dispersion could be assumed to be linear. For each of 
these subbands, a timeseries for all DMs up to the diagonal DM 
(i.e. the DM at which the dispersion delay across a frequency 
channel is about twice the sampling time) was formed. This re- 
sulted in a data set that can be thought of as a DM-frequency- 
time cube. 

At this stage the eight fan beams were combined into sub- 
beams by applying the appropriate geometric weighting over 
the total observing time before summing them, while still main- 
taining the DM-frequency-time cube structure. The appropri- 
ate weighting was determined by correcting for the fact that 
a particular position on the sky can move through different 
fan beams during an observation. The processing proceeded to 
cover the full range of the DM-frequency-time cube as follows; 
the dedispersion step was completed for 488 DMs separated by 
a DM step size of 0.696 up to the diagonal DM (of 340), result- 
ing in a timeseries of 2^^ samples for each DM. This time series 
was then Fourier transformed using a fast Fourier transform 
(EFT) routine. The resultant Fourier spectrum was cleaned by 
deleting frequencies known to correspond to interference and 
interpolated to recover spectral features lying between Fourier 
bins using the procedure described in [Ransom et al. (20021. 
This process was repeated for spectra in which 2, 4, 8 and 16 
harmonics were summed ( [Manchester et al.|[1985| l and for all 
DMs. All spectral features detected with a signal-to-noise ra- 
tio (S/N) of more than 7 were compared; those detected at mul- 
tiple DMs and apparently related harmonics were grouped be- 
fore selecting the combination of DM and spin frequency with 
the highest S/N for each candidate. 

For each of the resultant candidates, the partially dedis- 
persed data were used as the start of a refinement procedure 
where the data were dedispersed and folded at DM and period 
values around those found in the spectral search. The range of 
period and DM searched were determined by assuming a max- 
imum error of one Fourier bin. The partially dedispersed data 
were then divided into 8 to 32 subbands and 16 to 128 subinte- 
grations, depending on the spin frequency and DM of the candi- 
date, then folded and dedispersed at the period and DM values 
that gave the highest folded S/N, and stored for later inspection. 

The lists of potential candidates in all subbeams of one 
pointing were collated, and those with high S/N and detections 
in more than a single subbeam were stored. These detections 
were then compared to those made in the corresponding follow- 
up observation, and diagnostic plots, based on the aforemen- 
tioned analysis of the matches were inspected. The best candi- 
dates were reobserved in a follow-up timing programme, and 
three of those turned out to be pulsars. 



2.3. Survey Processing 

Apart from the initial processing required to form the eight 
individual subbeams, the survey processing used the same 
software as was developed for the Parkes multibeam surveys 



2.4. Follow-up timing 

The first follow-up series for the best pulsar candidates started 
in January 2006, and was again done in 8gr8 mode to be able 
to refine the positions. This also gave the advantage of cover- 
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ing a larger field-of-view, providing the possibility to use these 
observations as separate confirmation observations for any new 
candidate found in reprocessing of the original data. Sessions 
were separated by gaps, increasing in duration, to be able to 
phase-connect the calculated times of arrival for each obser- 
vation: in the first month of follow-up timing, observations 
were planned to take place in a sequence using day numbers 
{1,2,5,9,17,26). Then for half a year, monthly observations in 
8gr8mode were done, during which the positions were deter- 
mined with a precision accurate enough to be included in the 
normal WSRT timing programme. 



2.5. Monitoring sources 

When the position for each of the three sources was believed 
to be known with high enough accuracy, they were included 
in regular timing observations at the WSRT, using the nor- 
mal tied-array mode to optimise sensitivity. Furthermore, ob- 
servations were done at additional frequencies to determine 
the DM more accurately, and to determine the spectral indices 
of the pulsars. As part of the normal timing programme, the 
new pulsars were observed approximately monthly using PuMa 
( jVoute et all2 002). For observations at 328/374 MHz we used 
two lOMHz bands; the observations centred at 1380MHz or 
2270 MHz used 80 MHz of bandwidth, spread in 8 steps of 
lOMHz over a range of 160 MHz. The data were dedispersed 
and folded off-line, then integrated over frequency and time 
over the whole observation duration to obtain a single profile 
for each observation. Each profile was cross-correlated with 
a standard profile obtained via the summation of high signal- 
to-noise (S/N) profiles with the corresponding observing fre- 
quency (Fig.|5]), to calculate a time of arrival (TOA) for each 
observation. These were referred to local time generated by a 
H-maser clock at WSRT. The TOAs were converted to UTC us- 
ing global positioning system (GPS) maser offset values mea- 
sured at the observatory, and GPS to UTC corrections from the 
Bureau International des Poids et Mesures (BIPM|^ Finally, 
the TOAs were converted to the Solar system barycentre using 
the JPL epheme ris DE40f[^ We use d the new timing software 
package tempo2 ( Hobbs et al.|20061 l to analyse the data. 



3. Results 

The first three 8gr8 pulsars were discovered as relatively bright 
sources in the original 8gr8 observations. For each source, three 
observations were available at two or three different hour an- 
gles. As explained in Sect.|2] this allowed for removing any 
ambiguity in the position of the source, and measuring it accu- 
rately enough to observe the source in normal tied-array mode. 
To be able to obtain a coherent timing solution, a sequence of 



follow-up observations were done, as described in Sect. 2.5 



' http://www.bipm.org 

^ ftp://ssd.jpl.nasa.gov/pub/eph/export/DE405/de405iom.ps 
' http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Fig. 4. Period-period derivative diagram for the three discov- 
ered pulsars. The pulsars labeled with a square ar e confirmed 
pulsar-SNR associations (ATNF pulsar catalogu^ [Manchester 
let al.|2005D . 

3.1. PSRJ1937+2950 

PSR J1937-I-2950 is a 1.657 s pulsar with a characteristic age 
of about 8 Myr It was detected with high S/N in the original 
pointing and two confirmation observations. It was included 
in the first session of follow-up observations, and being de- 
tected in every observation, easily included in normal timing 
mode observations. The spectral index (SI) was measured to be 
-3.3(2), calculated using the procedure described in [Lorimer 
et al. ( 1995| l. The pulsar was not detected at 2.3 GHz, which is 
not surprising given the large SI. Two observations at 140 MHz 
failed to detect the pulsar despite its very large spectral in- 
dex. Interstellar scattering effects for this pulsar are expected 
to be negligible ( Cordes & Lazio|2002 i. The non-detection at 
140 MHz may therefore indicate that the spectrum has flattened 
or turned over in this range. Summed profiles for the observed 
frequencies are shown in Fig. [5] and the timing solution is pre- 
sented in Table[T] 



3.2. PSR J2027+4557 

PSR J2027+4557 has a period of 1.099 seconds. It is relatively 
old (57 Myr). It was independently detected in the GBT350 sur- 
vey at Green Bank (Hes sels et al.||20(j8] l. The timing solution 
is presented in Table[T| The pulse profile exhibits quite a bit 
of scattering at lower frequencies, see Fig.|5] which is consis- 
tent with that expected for the line of sight towards this pulsar 
( [Cordes & Lazio|2002| l. The SI of this pulsar was measured to 
be -1.43(9), based on measurements of the mean flux density 
from observations at 328, 374, 1380, and 2300 MHz. 

Of the three pulsars we discovered so far, PSR J2027+4557 
is the brightest. Its individual radio pulses are easily detected, 
allowing for an analysis of subpulse modulation. The shape and 
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Fig,5. WSRT standard profiles for PSR J1937+2950 (left), PSR J2027+4557 (middle) and PSR J2047+5029 (right). These tem- 
plates have been generated using ~ 15 hours of timing data, except for the profiles of PSRs J2027+4557 and J2047+5029 at 
2300 MHz, where ~2 hours of timing data was used. 



intensity of the individual pulses is variable from pulse to pulse, 
and this can be visualised by making a pulse-stack: a plot in 
which successive pulses are plotted on top of each other. The 
modulation patterns in such plots can be extremely irregular, 
or very regular, depending on the pulsar. Repeating patterns 
of diagonal bands of emission, the so-called drifting subpulses 
( Drake & Craft 1968[), are observed for many pulsars in varying 



degrees of clarity (Weltevrede et al. 2006 Weltevrede et al 



2007| l. Subpulse modulation can be analysed by calculating the 
longitude resolved fluctuation spectrum (LRFS; [Backer|1970)l 
and the two-dimensional fluctuation spectrum (2DFS; 'Edwards' 



& Stappers 2002 1. For a detailed description of the methods 
used below, we refer to |Weltevrede et al.] ( |2006| l. 

Visual inspection of the individual pulses of 
PSRJ2027+4557 reveal that pulses at the leading edge 
become stronger and weaker on a timescale of slightly more 
than 10 pulse periods, a clear sign that this pulsar could have 
drifting subpulses. The result of our fluctuation analysis of 
this pulsar can be seen in Fig. [6] The top panel shows the 
pulse profile, and the longitude-resolved modulation index is 
superimposed. The profile shows at least three components (at 
158°, 161°and 172°pulse longitude). The LRFS (middle panel) 
shows, for each pulse longitude, the Fourier transform of the 
intensities of the successive pulses, allowing us to determine 
whether the modulation is random or periodic. The most 
prominent feature is the patch of power associated with the 
leading component with a quasi-periodic frequency of about 
0.07 cycles per period, thereby confirming the aforementioned 
intensity modulation with a period P3 = 14 + 1 pulse periods. 
The lower panel shows the 2DFS of the leading component, 
which is sensitive to the slope of the drifting subpulses. The 
centroid of the power in between the two dashed lines is 



offset from the vertical axis, indicating that the subpulses drift 
towards later times in the pulse-stack. Other parts of the profile 
do not show evidence for drifting subpulses. 

3.3. PSRJ2047+5029 

PSR J2047 -1-5029 was first detected at its second harmonic. 
Separated by approximately 180 degrees, the pulsar has an 
interpulse that is about 50% as strong as the main pulse. 
Like PSRJ2027-H4557, PSR J2047H-5029 is relatively bright at 
21 cm, having an SI of -1.39(9). It shows quite a bit of timing 
noise, see Fig.|7] There is a weak hint of periodicity in the tim- 
ing noise, however it is too early to determine if there is any 
physical intrinsic cause for the timing noise, like a long-period 
binary orbit or small glitches. 

Association witli supernova remnant HB21? 

PSR J2047 -1-5029 is located along the same line of sight as the 



supernova remnant (SNR) HB21 (Hanbury Brown & Hazard 
[1953, G89.0-I-4.7), with its position ~ 0.5 degrees away from 
the centre of the SNR as given by Green's catalogu^ see 
Fig.|8] HE 21 is irregularly shaped, probably due to a molec- 
ular cloud on one side that may have stopped the remnant from 
expanding symmetrically ( Tatematsu et al.|1990 1. 

HB21 was searched for a radio pulsar before (Biggs & 
Lyne|1996| , but down to a Hmit of 13 mJy (at 610 MHz) nofli 



ing was found. The search covered an area 0.5 degrees around 
the SNR centre, which probably just included the pulsar po- 
sition. However, interpolating from our SI calculations, the ex- 
pected flux density at 610 MHz is 1.23 mJy, explaining why the 

■* http://www.mrao.cam.ac.uk/surveys/snrs/ 
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Pulsar name 


J1937+2950 


J2027+4557 


J2047+5029 


Fit and data-set 


MJD range 


53463—54751 


53462—54751 


53740—54751 


Number of TO As 


85 


89 


101 


Rms timing residual (ms) 


2.8 


0.6 


1.3 


Weighted fit 


Y 


Y 


Y 


Reduced value 


1.2 


1.7 


5.5 



Measured Quantities 



Right ascension, a (J2000) 

Declination, S (J2000) 

Pulse frequency, v (s"') 

First derivative of pulse frequency, v (s^^) . . 
Second derivative of pulse frequency, v (s^^) 
Dispersion measure, DM (cm^^^ pc) 



19*'37"47>;603(14) 

+29°50'0P:8(2) 
0.603344173844(9) 
-1.2672(6)x 10-15 

113.80(2) 



20''27""lffi233(3) 
+45°57'57'.'08(4) 
0.909379126021(5) 
-2.559(3)xlO-'« 

229.594(11) 



20''47"' 54^.6400(4) 
+50°29'38'.'17(4) 
2.242431224778(16) 
-2.09976(8)xl0-''i 
-7.41(8)xlO-2-+ 
107.676(5) 



Set Quantities 


Epoch of frequency determination (MJD) . . 
Epoch of position determination (MJD) .... 
Epoch of DM determination (MJD) 




54250.0 
54250.0 
54250.0 






Derived Quantities 


Characteristic age, t (Myr) 


7.5 




56.3 


1.7 


Surface magnetic field strength, S (10'- G) 


2.4 




0.6 


1.4 


Distance, d (kpc) 


5.1(7) 


10.5 (7.4-50.0) 


4.4(5) 


Energy loss, (erg s"') 


3.0x10^' 




9.2x10'" 


1.9x10" 


Energy loss/d-, E^aJd- (ergs"' kpc"^) 


1.2x10^° 




8.3x102** 


9.6x10^' 


Mean flux density at 328 MHz, (mJy) 


8.7(6) 




9.1(9) 


2.5(3) 


Mean flux density at 1380 MHz, (mjy) 


0.08(2) 




1.34(13) 


0.38(4) 


Mean flux density at 2300 MHz, (mJy) 






0.32(9) 


0.22(13) 


Spectral index 


-3.3(2) 




-1.43(9) 


-1.39(9) 


Assumptions 


Clock correction procedure 

Solar system ephemeris model 

Model version number 




TT(TAI) 
DE405 
5.00 







Note: Figures in parentheses are the nominal Icr tempo2 uncertainties in the least-significant digits quoted. 
The DM distances are estimated from the |Cordes & Laziol ( |2002^ model. 



pulsar was not detected in their survey. 'Lorimer et al. ( 1998) 1 
searched the whole remnant, also at 6 10 MHz, to a sensitivity 
limit of 0.66 mJy (this number is calculated accounting for the 
duty cycle of the pulsar of ~0.02). Somewhat surprisingly, they 
did not detect PSR J2047-H5029. 



In general, the probability of having a pulsar along the line 
of sight towards an SNR by coincidence appears to be high, 
considering the number of pairs that are shown to be not associ- 
ated (e.g. jGaensler & Johnston] 1995) . Recent targeted searches 
have revealed a handful of pulsars in SNRs, but only a minor- 
ity of those are considered to be real associations (e.g. [Camilo 
et al.||2002[ [Lorimer et"aL]|2006[ ). |Kaspi| ( |1996 ) has proposed 
a number of questions to assess whether a pulsar is associ- 
ated with an SNR, which we will address one by one below 
for this PSR-SNR combination. From the following, we con- 
clude it is very unlikely that PSR J2047+5029 is associated 
with SNR HB 21; 



Distance: Estimates of the distance of SNR HB21 vary a lot. 
References in |Leahy| ( |1987) l give a range of 1.0 to 1.6 kpc, all 
derived using the surface biightness - diameter (S - D) rela- 
tion for SNRs. From radio continuum, HI and CO observations. 



Tatematsu et al.| \\99Q) state that the SNR is part of the Cyg 
OB7 association, which has a known distance of 800(70) pc 
( [Huang & Thaddeus|1986 i However, due to the lack of direct 
evidence of interaction of the SNR with the OB association, 
combined with a low X-ray surface brightness, the distance to 
HB21 was proposed to be larger than 1.6 kpc (references in 
[Byun et al.|2006| l, and they combined all available distance es- 
timates to arrive at a value of 1.7(5) kpc. 

From the measurements of the DM of PSR J2047+5029, 
using the jCordes & Lazi"o ( 2002| l model, we derive a distance 
of 4.4(5) kpc. Although this may be overestimated due to un- 
modelled features in the DM model, the distance estimates for 
pulsar and remnant are at least different by a factor 3. However, 



in Fig. 6 of Cordes & Lazio (20021, a schematic overview is 
given of location of spiral arms in the Galaxy, which suggests 
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1 40 1 60 1 80 

Pulse longitude (deg) 
I I I I I .. I I I I I I I I 




0.02 



-100 100 

Fluctuation frequency CPq/P2) 

Fig. 6. Top panel: Pulse profile and Longitude resolved mod- 
ulation index for PSR J2027+4557. Middle panel: longitude 
resolved fluctuation spectrum (LRFS). Bottom panel: two- 
dimensional fluctuation spectrum (2DFS). These plots show 
one 50min. observation made using the new pulsar machine, 
PuMall (Karuppusam y et al.|2008| l, which provided twice the 
bandwidth at 1380 MHz and thereby increased the sensitivity 
to single pulses. 



that the pulsar, if at at distance of 4.4 kpc, is located in (or be- 
hind) the Perseus arm. This argues in favour of a coincident 
association. 




2006 



2007 



2008 



2009 



Year 



Age: The age of the SNR was estimated by |Leahy| ( |1987| l at 
8000-15000 yr, by comparing the properties of the SNR (shock 
temperature, interior density) with models for SNR expan- 
sion into a 3 -component medium. Leahy & Aschenbach ( 1996) 
quote an age of 19000 years, based on a model for an SNR in a 
cloudy interstellar medium. 



Fig. 7. Timing residuals for PSR J2047-H5029. Top: solution in- 
cluding astrometric parameters, spin period and first spin pe- 
riod derivative. Bottom: Solution also includes second deriva- 
tive. Different symbols represent different observing frequen- 
cies; circles: 323 and 374 MHz, squares: 1380 MHz. For the 
parameters of this solution, see Table[T] 

The characteristic age of the pulsar is two orders of mag- 
nitude larger, about 1 .7 Myr. This apparent age conflict is yet 
another argument against association of the pulsar and the rem- 
nant. The characteristic age of a pulsar must be regarded as 
effectively only an upper limit to its real age, when a dipole 
braking model is assumed. However, requiring the pulsar age 
to be consistent with that of the SNR age and assuming the 
pulsar period derivative has not changed during its lifetime, 
leads to the conclusion that the biith period of the pulsar would 
have been similar to the cuiTent period to the extent of 10"^. 
As almost all pulsars that are confirmed to be associated with 
an SNR have spin periods less than about 140 ms ( [Lorimer 
|et al. 11998) [Kramer et al.|2003] l, the measured spin period for 
PSRJ2047-H5029 of 0.445 s is very unlikely to be close to its 
spin period at birth. So, though the characteristic age may not 
be used as a realistic age estimate, the pulsar must be quite 
old to have slowed down to this period, and we conclude that 
PSR J2047-H5029 is too old to be associated with the SNR. 

Velocity/Proper Motion: Most pulsars receive a kick velocity 
when born in a supernova event. The angular distance between 
the current position of a pulsar and the centre of its claimed 
corresponding SNR, combined with the estimated age and dis- 
tance towards the combination, can therefore be used to calcu- 
late a minimum velocity needed for the pulsar to travel to its 
present position in its lifetime (e.g. Migliazzo et al.|20"02 l. 
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neous signal-to-noise squared with the field-of-view, or, 



20"48'" 20"44'" 
Right Ascension (J2000) 



Fig. 8. hnage of SNR HB 21, using WENSS radio images 
(jRengelink et al.|[T997] l, and ROSAT contours ( |Leahy 
|Aschenbach"1996 



& 



The symbols represent various centre po- 
sitions of the SNR as quoted in the literatu re. (-H,Green,,2006] 
Leahy & Aschenbach||l996| and x |Leahy||1987| The ★ rep- 



resents the best position of PSRJ2047 -1-5029 from timing, see 
Table[T] The error on pulsar position is smaller than the size of 
the symbol. 



For PSR J2047+5029, depending on the estimated age and 
distance of the remnant, the required velocity lies between 300 
and 2000km s"', which is large but not unreasonable. So far we 
have not measured a proper motion signal in the timing resid- 
uals. Given the timing precision that we can achieve for this 
pulsar, and as we have been timing the pulsar for only a lim- 
ited timespan, this is not surprising because the expected effect 
that a transverse velocity of 300-2000 km s ' would have on the 
timing residuals is not measurable. It is therefore at this stage 
not possible to use proper motion as an argument for or against 
association between HB 21 and PSR J2047+5029. 



Interaction with tine remnant: There is no indication of emis- 
sion at the pulsar position that might be from the pulsar it- 
self, for a pulsar wind nebula, or from interaction with the 
SNR itself (See Fig.[8]|. However that may be expected from 
the flux limits of the WENSS survey (Renge Unk et al.|1997" i. 
Unfortunately there is no radio image available with sensitiv- 
ity to have detected the pulsar. Therefore this criterium is not 
useful for determining whether the pulsar and remnant are as- 
sociated. We note also that there is no indication of either a 
pulsar or an associated wind nebula in the X-ray emission. 



4. Pulsar Surveys with Interferometers 

As previously discussed in ^ interferometers have not gen- 
erally been used for pulsar surveys due to their limited survey 
speed, SpsR, which can be expressed as the product of instanta- 



2PSR = (^) v^'^WQ,., 



sys 

where Aeff is the effective collecting area in m^, Tsys the system 
temperature, v the observing frequency in GHz, a the typical 
pulsar spectral index of about -1.8, BW the observing band- 
width in GHz and Q^z the field-of-view in square degrees. In 
view of the steep spectral indices of many pulsars, there is a 
natural speed advantage for moderately low frequency surveys, 
but there are often other considerations which dictate the opti- 
mum choice. In any case, this metric can be used to compare a 
number of existing and potential future pulsar surveys as shown 
in Table|2j where we note that higher values represent higher 
survey speeds. 

The first two entries in the table document the very suc- 
cessful pulsar surveys carried out with the Parkes telescope at 
70 and 20 cm ( Manchester et al.|[T996l |2001| l. A similar high 
survey speed was achieved in both bands by employing both 
multiple beams and much larger bandwidths at 20 cm relative 
to 70 cm. The current Arecibo 20 cm multi-beam ALFA system 
( [Cordes et al.|2006) l offers a substantial improvement in survey 
speed over these previous efforts. In contrast, the current WSRT 
8gr8 survey already achieves about twice the ALFA pulsar sur- 
vey speed while employing only 10 MHz of bandwidth (while 
20 MHz can now be accomodated with the Pumall upgrade in 
place). 

The final three lines in the table offer guidelines for what 
will be necessary for future interferometric systems if they are 
to provide competitive pulsar survey capabilities. The ASKAP 
(Australian Square Kilometre Array Pathfinder) telescope is 
now under construction on the short-listed SKA site in Western 
Australia ( (Johnston et al.|[2007| l. It will consists of 36 tele- 
scopes of 12m diameter, fed in the 0.7 - 2 GHz band by a fo- 
cal plane array feed that will provide good performance over a 
30 square degrees field-of-view. The ASKAP configuration is 
somewhat centrally concentrated but is distributed over a re- 
gion of about 6 km in extent. In the table we have attributed 
only 50% of the total ASKAP sensitivity to a synthsized beam 
of about 100" diameter at 850MHz. Since the ASKAP con- 
figuration is a randomised two dimensional one, rather than a 
pure "grating array" such as the WSRT, individual tied-array 
beams will need to be generated for each search direction in 
real time. Competitive pulsar surveys with ASKAP will require 
something like 8000 of these simultaneous tied-array beams to 
provide a significant improvement in survey speed. Together 
these would cover only about 6 of the 30 square degrees total 
field-of-view that is available. The meerKAT project plans to 
construct an SKA pathfinder in the other short-listed SKA site 
in Southern AfriccPl It will consist of 80 telescopes of 12 m di- 
ameter, fed in the 0.75 - 2.3 GHz band by a single pixel feed. 
The configuration is such that 70% of the collecting area will be 
contained in the inner 700 m. As with ASKAP the non-grating 
array layout will require beam formation to be done in real 
time. The larger collecting area concentrated into a smaller re- 
gion mean that meerKAT would require fewer beams, about 

http ://www. ska. ac . za/meerkat/index . shtml 



10 



G. H. Janssen et al.: Discovery and timing of the first 8gr8 Cygnus survey pulsars 



Table 2. Pulsar survey speed comparison. 



Pulsar Survey 




V 


BW 




SpSR 




(m^/K) 


(GHz) (GHz) (deg2) 




Parkes 20 cm Multi-beam 


70 


1.4 


0.30 


0.41 


180 


Parkes 70 cm 


32 


0.41 


0.03 


0.39 


290 


Arecibo 20 cm ALFA . . . 


805 


1.4 


0.3 


0.02 


1100 


WSRT90cm 8gr8 


31 


0.33 


0.01 


4.2 


2230 


ASKAP 35 cm 8k-beam . 


32 


0.85 


0.3 


6 


3300 


meerKAT 30 cm Ik-beam 


180 


1.0 


0.5 


0.5 


8300 


SKA 35 cm Ik-beam .... 


2000 


0.85 


0.3 


0.75 


1.6x10* 



400, to be a competitive search instrument. For Table|2] we 
give numbers for a 1000 beam survey (J.Jones, priv.comm.). 
The much higher sensitivity expected for the SKA, of about 
2000m"^/K within the central few km, permits very compet- 
itive pulsar surveys with many fewer simultaneous tied-array 
beams. For example, 1000 beams would suffice to cover about 
0.75 square degrees and yield about 3 orders of magnitude im- 
provement in survey speed over the current state-of-the-art. We 



direct the reader to Smits et al. ( 2009 1 for a more complete dis 



cussion of SKA surveys. 
5. Conclusions 

Our new pulsar survey technique at the WSRT has resulted in 
the discovery of three new pulsars. In the first analysis of the 
data three new pulsars have been found. In contrast to the ex- 
pectation of finding young pulsars, these are normal pulsars. 
However, the analysis procedure is being optimised and we ex- 
pect finding many more pulsars and transients in the Cygnus 
region. The results and discussion of the full analysis of the 
survey will be presented in a forthcoming paper. 

One of our newly discovered pulsars, PSR J1937+2950 has 
an SI of -3.3(2), showing that the survey is successful in find- 
ing steep-spectrum pulsars, and that there are potentially more 
to find. Somewhat surprisingly, two of the pulsars are brighter 
at 2 1 cm, indicating that surveys of this region at that frequency 
may also still find quite bright pulsars. Moreover, the discov- 
ery of PSR J2027 -1-4557 shows that it is still possible to find in- 
teresting single pulse sources with modern surveys. Following 
several evaluation criteria, we conclude that PSR J2047+5029, 
although having a position along the same line of sight as 
SNR HB 21, is not associated with this SNR. 
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